An exactly solvable dead- layer problem 
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A molecule consisting of two particles interacting with a delta-function attraction, and confined to a one- 
dimensional volume, is studied. From the exact solution of the system we deduce that the center-of-mass 
effectively moves in a volume reduced by an inaccessible dead layer near each wall. For a very large volume 
the dead-layer width equals the size of the molecule, defined as the expectation value of the interparticle 
distance for the unconfined molecule. For finite volumes the dead-layer width is smaller than this. For 
unequal particle masses, perturbation in the mass difference up to third order yields the same result, and 
we conjecture that for this system the dead-layer width is independent of the mass values. 

Keywords: Two-particle system, exact solution, dead layer. 

1 Introduction 

For the center of a small classical ball in a container, there will be an inaccessible volume near 
the wall. Should the ball be incompressible, the width of the inaccessible volume is equal to the 
ball radius. A similar situation is expected for a bound quantum-mechanical system, a molecule 
or an exciton, say, in a large finite volume. The center of mass cannot easily reach the boundary 
since the separate particles have to stay inside the volume. Therefore the center of mass effectively 
moves in a volume that compared with the actual volume is reduced by a dead layer at the bound- 
ary. Excitons or other bound systems in quantum dots or quantum wells are prime examples (Fig. 1). 

Consequently one expects the ground-state energy of the confined system to be of the form 

h 2 ir 2 



for a large spherical volume of radius R, and 



fc2 2 

n ir 



E =- B+ 2M(L-^ (2) 

for a quantum well of width L. Here B is the binding energy of the system when unconfined, M 
is the total mass of the system, and d measures the width of the dead layer. The notion of a dead 
layer was apparently first used by Pekar [TJ. The energy expressions ((T|) or (|2j) arc commonly used 
by experimentalists to fit their data [H [3J 0] . 

From a fundamental point of view the dead-layer concept raises several questions: 

• A proper dead layer should be independent of the volume in the limit of a very large volume. 
While one may always write the energy in the form ([TJ or (|2|) , one should show that the value 
of the dead layer parameter d approaches a limiting value when the volume increases. 



* Corresponding author. Email: chr-hem@online.no 



1 



(a) 



(b) 



Figure 1: Dead layer in a spherical quantum dot of radius R, and in a quantum well of width L. 



• It would be interesting to see whether d approaches its limiting value from above or from 
below. Some experimental results [3] indicate that d is smaller than its limiting value when 
the volume is finite. 

• The size of an unconfined system of two particles with masses mi and 7712 is determined by the 
reduced mass /i = 77117712/(7711 + 7712). However, there is no reason to expect that the dead-layer 
width depends on the masses merely through fj,. Whether, for a given value of the reduced 
mass, the width of the dead layer increases or decreases with increasing difference between the 
two masses is controversial (see [4] and references therein). 

We will study these questions in a confined two-particle model. In order to have an exactly 
solvable system we work in one dimension, and with an attractive delta-function interaction between 
the particles. For equal masses we will be able to find the exact ground state of the confined system. 
This is of interest in itself as a new non-trivial solution of the Schrodinger equation. In Sec. 2 this 
exact solution is worked out. The effect of nonequal masses is studied in Sec. 3, and our results are 
summerized in Sec. 4. 



2 The model and its exact solution 

We consider two particles with delta-function attraction in a one-dimensional quantum well of ex- 
tension L with infinite offsets. The Schrodinger equation is 

n 2 d 2 i> n 2 d 2 i, h 2 

6{X2 — Xx)tp = Eyj. (3) 



2mi dx\ 2ui2 dx\ fia 



Here fi = TOi?7i2/(rni + TO2) is the reduced mass, and a is a measure of the potential strength. The 
particles live on the interval (— L/2, L/2). Hence the boundary conditions imposed on the wave 
function ip(xi, X2) by the rigid walls are 

^(xi,±L/2) = V(±L/2,i a ) = 0. (4) 

The Schrodinger equation ([3]) can alternatively be expressed in terms of the relative coordinate 
x = x-x — X\ and the center-of-mass coordinate X = (rn\Xi + 1712X2) j (mj + m-i) as 

h 2 d 2 ^b h 2 d 2 %b h 2 

V V 6(x) ^ = £V, (5) 



2M dX 2 2n dx 2 /ia 

with the total mass M = mi + m 2 . 

In the absence of walls the interaction between the two particles allows the bound state 

4>(x) = a~ 1/2 e~ M/a , (6) 

with binding energy 

The mean linear size of the bound state is of order a, 

(\x\)=a/2. (8) 

The rms size is a factor \[2 larger. 

The important parameter of the problem is the dimensionless ratio a/L. For the dead-layer 
problem we are interested in situations in which the molecular size a is much smaller than the 
container size L. 



2.1 The exact solution for equal masses 



We now assume mi — m\ — m. Since the Schrodinger equation is then symmetric in the position 
variables x\ and xi-, the ground state must also possess this symmetry. 

We first solve the Schrodinger equation in the two domains x\ < xi and x\ > X2, in which the 
pair potential does not appear explicitly. The two solutions are then matched to ensure that the 
wave function along the diagonal x\ = X2 is continuous and has the singularity appropriate for the 
delta-function interaction. 

In the domain a; 1 > X2, the appropriate boundary conditions are ip(L/2, x^) = and ip(xi, —L/2) = 
0. The Schrodinger equation ([3]) may be solved by separation of variables. Product solutions satis- 
fying the boundary conditions are 



^1^2(^1 1^2) oc sinh[Ai(a;i - L/2)/L] sinh[A 2 (i£2 + L/2)/L]. 
The two dimensionless separation constants Ai and A2 satisfy 



(9) 



E 



(Xj + X 2 ). 



2mL 2 

By symmetry we have in the other domain x± < X2 ■ 

i>\ u \ 2 (x\,X2) oc sinh[Ai(x 2 - L/2)/L] sinh[A 2 (a;i + L/2)/L]. 



(10) 



(11) 
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By construction the wave function ip\ lt \ 2 ^ s now everywhere continuous. We may write more com- 
pactly 

"0Ai : a 2 <* smh[X 1 (x 1 +x 2 + \x 2 -x 1 \- L)/2L] sinh[X 2 (x 1 +x 2 -\x 2 -x 1 \ + L)/2L] (12) 
= sinh[Ai(2X+ |x| -l)/2] sinh[A 2 (2X - |x| + l)/2] (13) 
oc - cosh[(Ai + \ 2 )X + (Ai - A 2 )(|a| - l)/2] 

+ cosh[(A 1 - \ 2 )X + X 1 + X 2 )(\x\ - l)/2]. (14) 

We have introduced dimensionless relative and center-of-mass coordinates, 

X = (zci + x 2 )/2L and x = (x 2 — xi)/L, (15) 

and used that 2 sinh(u) sinh(u) = cosh(w + v) — cosh(u — v). 
The superposition 

iP(X,x) = Ai> XlM (X,x) + Bip X2tXl (X,x), (16) 
with arbitrary constants A and B, is also a solution, since both terms correspond to the same energy 

It is clear that the separation parameters Aj cannot always be real, since for sufficiently weak 
interaction the energy must be positive, corresponding to imaginary Aj. We therefore take Ai and X 2 
complex. Since the energy (jlOp is real, we select separation parameters that are complex conjugates. 
Thus we put 

Ai = u + iv and A2 = u — iv, (17) 



with real u and v. Moreover it is clear that with A = B in (|16|) the wave function will be real, and 
we will show that this is the correct superposition for the ground state. Using (|14|) . (fT7| and (|16|) 
with A — B, we obtain 

ip(X,x) oc ^Re[cosh(2ivX - u + u\x\) - cosh(2uX - iv + iv\x\)] (18) 
= cos(2wX) cosh(u|x| — u) — cosh(2uX) cos(ii|x| — v) (19) 

By construction this wave function obeys the boundary conditions, and also the Schrodinger 
equation for x\ ^ x 2 . It remains, however, to show that it may satisfy the Schrodinger equation 
with the correct interaction potential V(x) = —5(x)h 2 / /j,a. 

For a general interaction potential V the Schrodinger equation in dimensionless relative and 
center-of-mass coordinates takes the form 

U 2 d 2 ip H 2 d 2 ip , N 

-2M^^-W^ + ^ = ^- (20) 

Thus the potential can be calculated from the wave function, 

n 2 ay n 2 ay 

2fiL 2 V dx 2 2ML 2 ijj dX 2 1 ' 

Clearly the desired delta-function potential can only come from the first term containing a differen- 
tiation with respect to the relative coordinate x. Using d 2 \x\/dx 2 = 28{x) we obtain 

h 2 v sin(u) cosh(2itX) + u sinh(w) cos(2uX) 
~~ ^2 cos(u) cosh(2uX) - cosh(u) cos(2t>X) ^ ' 
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The desired form in ([3]), 



h 2 , h 2 



V= 5{x 2 -x 1 ) = -5{x), (23) 

jia fiaL 



with dimensionless x, will be obtained when we take 

usin(w) = cos(v) (24) 

a 

and 

usinh(u) = — cosh(w). (25) 
a 

In conclusion, we have found the following exact ground state of the system: 

ip(X, x) = iV[cosh(w|a;| - u) cos(2vX) - cosh(2uX) cos{v\x\ - v)], (26) 
where the coordinates are scaled by L. The parameters u and v satisfy 

utanh(u) = — vtan(i;) = L/a, (27) 

and the energy (ITU1) is given by 

E=-^(v 2 -u 2 ). (28) 
It is straightforward, but tedious, to determine the normalization constant N. The result is 




(29) 



(sinh(2u) + 2u)(sin(2u) + 2v) 



2.2 Discussion of the solution 

Both the parameters u and v may be taken nonnegative, since the wave function is invariant under 
parameter sign changes. Since u tanh(u) is a monotonically increasing function, u is determined 
uniquely by ([27j) . The equation for v, on the other hand, has no unique solution. Since the energy 
expression ([28} increases with increasing v, the ground state corresponds to the smallest solution 
for v. (|2"T|) shows that tan(u) is negative, and we conclude that the relevant solution for v lies in the 
interval (ir/2,n). 

The boundaries of the w-interval correspond to special limiting cases. For v = it we have u — 
and 

2 

tp = L [cos(2-7rX) + cos(7r|x|)] = — cos(nxi/L) cos(ttx2/L) 7 (30) 

±j 

two uncoupled particles in the domain. This strong confinement limit is consistent both with the 
corresponding energy value, and with (|27|) , since v — ir corresponds to a = oo, i.e., vanishing 
delta-function amplitude. 

When v — > tt/2 from above, eq. (f27|) shows that u becomes very large approaching the value 
L/a. This is the weak confinement limit, and may be considered as a being finite and L — > oo. The 
interesting factor in the wave function (|26[) is 

cosh(u|a;| — it) ~ cosh(u) cosh(u|x|) — sinh(it) sinh(it|x|) ~ cosh(it)[cosh(it|o;|) — sinh(u|a;|)] 

= cosh(it) e~ ulxl ~ cosh(u)e- lx2 - x ^ /a . (31) 
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Using also that the normalization factor is approximately given by 1/N ~ J aL/2 cosh(u) in this 
limiting case, we obtain altogether 

^~a- 1/2 e- |a;2 ^ l|/a (2/L)- 1 / 2 cos{ir(x 1 +x 2 )/L), (32) 

a product of the bound state and the center-of-mass ground state in a large volume. 

To illustrate the wave function inbetween these limiting cases we show in Fig. 2(a) tjj(x,0) for 
different ratios a/L. The cusp at x = becomes more and more pronounced with increasing L/a, as 
expected. The second derivative near the cusp is positive only for a < 0.431L. In Fig. 2(b) tp( x i> ^2) 
is shown for L — 5a. We note that the wave function is everywhere non-negative, as must be in the 
ground state. 



2.3 The dead layer 

We now investigate the energy (l28l) for large volumes. For L/a large we see from (|27p that u is large 
and approximately equal to L/a. We take u = L/a, since the correction terms to this value are 
exponentially small and can safely be ignored. Thus we have to discuss 

By eq. the first term is the molecular binding energy —B. Writing the energy (|33|) in the 
suggestive form 

£?= - g+ 2M(L-2d )a ' (34) 

we have 

L ttL 
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Figure 3: The effective dead iayer as function of L/a. 

Here v € (ir/2,ir) is given by wtan(w) = —L/a, eq. (|27|) . It is straightforward to solve this equation 
for v perturbatively for small a/L, since v is very close to tt/2 in this limit. We find 



7T 

V= 2 



a a 2 ( tt 2 \ a 3 / a 4 



Insertion into ea. (|35|) gives the corresponding expansion for d: 

a 7r 2 a 3 

Thus in the limit of a large container a dead layer of width 

d = a/2 



(36) 

(37) 
(38) 

of the unconfincd 



exists. Hence the dead-layer width equals precisely the mean linear size 
molecule. 

Moreover, when the container is finite, the effective dead layer is somewhat reduced in size. This 
can be see in more in more detail by calculating d for all values of a/L. We must then use the full 
equation (|27[) for u, not the weak confinement approximation u ~ L/a. Fig. 3 shows that d is a 
monotonically increasing function of L/a. 

The molecule is expected to be deformed close to a wall. To illustrate this, we calulate the linear 
size of the molecule, 

J\x\ \ip(X,x)\ 2 dx 



<M> = 



f\i>{X,x)\ 2 dx ' 



(39) 
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Figure 4: The dimensionless molecular size (\x\}/a as function of (L/2 — X)/a, the dimensionless 
distance between the center-of-mass position and the boundary. 



for X positive, and study how the size changes when the center of mass approaches the wall at 
X = L/2. It is straightforward to show that for large L/a, the molecular size depends merely upon 
L/2 — X, the center-of-mass distance from the boundary, independent of the container size. The 
integration in (JSJ yields for large L/a: 

„ i \ tf z [(y-l)e- 2z -(2z-l)e-y} 2 ydy , 

fi z [(y-l)e-**-(2z-l)e-y}2dy 
a (2z - l) 2 - (16z - 8)e- 2z + (16z 4 + f z 3 + 20z 2 - 9)e~ 4z 



where 



2 (2z-l) 2 + {f z 3 + 4z 2 -l)e- 4z 



z=(L/2-X)/a. (42) 



Details of the integration are not given. For large z — > oo, i.e. far from the wall, we recover from 
(|4T|) the molecular size (|a;|) = a/2, as expected. Fig. 4 shows how the molecule is compressed when 
its center of mass is forced to be close to a wall. When the center of mass is at a dead-layer distance 
from the wall, X = L/2 — a/2 or z — 1/2, the molecular size has by eq. (|4T]) shrunk to one-half of 
its size in bulk. 
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3 Unequal masses 



For particles with unequal masses the hamiltonian is 



2m i dx\ 2m,2 dx\ afi^ X2 ^ 

The boundary conditions are as before. In the absence of an exact solution of the correspond- 
ing Schrbdinger equation, we now resort to Rayleigh-Schrodinger perturbation theory in the mass 
difference. 

3.1 First-order perturbation 

We use an unperturbed Hamiltonian 

h 2 d 2 h 2 d 2 h 2 

Ho = ~T t: 2 7\ "o 2 5{x 2 -xi), (44) 

2m ax\ 2m axr^ a/i 

with m selected such that the reduced mass /i is the same for the perturbed and the unperturbed 
problem: 

A--L + -L-I. (45) 

m mi m2 ji 



The perturbation is then 



H' = H-H ='-L(^-^)(^-^). (46) 



/ j_ j_\ _ d 2 

4 \m\ TO2 / \dx\ dx\ 

The ground state (|26p of the unperturbed problem (I44[) is symmetric with respect to interchanging x\ 
and X2, while the perturbation H' is unsymmetric. Hence the hrst-order perturbation contribution 
vanishes: 

= (0|ff'|0) =0. (47) 

3.2 Second-order perturbation 

The second-order perturbation contribution to the ground state is 

i^V'T'T (48) 
^ El -El 



To evaluate the sum we use the Dalgarno-Lewis technique [5]. Suppose we can find an operator F 
satisfying 

H'\0) = {FH - H o F)\0) = {F,H ]\0}. (49) 

Then 

(k\H'\0) = (E%-E%)(k\F\0), (50) 
and the sum in (|48j) may be written as 

Ei 2) = J2 1{ ^ H J°J = J2(0\H'\k)(k\F\0) = (0\H'F\0) (0\H'\0){0\F\0). (51) 

k^O k k^O 
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The great advantage is that the summation is avoided, and that merely the unperturbed ground 
state enters. 

In relative and center-of-mass coordinates, we have 



H 



d 2 



n 2 d 2 



m dx 2 4m dX 2 afi 



and 



H' 



1 



1 



2 V mi m 2 J dx dX 



d d 



For the problem at hand the Dalgarno-Lewis operator F may be taken as 



1 



1 



F=- [ 

mi m 2 



d 

dX' 



(52) 
(53) 

(54) 



In (|49|l the commutator between F and Hq enters. Since d/dX commutes with Hq, eq. (|52|) . we have 



[F,H }=-—[ 

mi m 2 



x, 



dx 2 



d 

dX 



h 2 



i 



l 



2 V mi vni J dx dX 



d d 



(55) 



thus the Dalgarno-Lewis requirement (|49|) is fulfilled [6]. 

It only remains to determine the right-hand side of eq. ([51]) 



n m 



mi m.2 



9 a 2 







(56) 



dx dX 2 

The evaluation of the ground state average (f5"6"f is tedious and is referred to the appendix. We obtain 



n 2 7r 2 


( m 


--) 2 




lQmL 2 


l mi 


m 2 / v 


+ 2 l) 



(57) 



to first order in the small parameter a/L. 
3.3 Third-order perturbation 

The third-order perturbation contribution to the ground-state energy would depend on the masses 
mi and m 2 as follows 

EM<x(- -V. (58) 

\mi m 2 / 

But it would be unphysical if the energy contribution would change sign by interchanging the mass 
values. Thus the proportionality constant must be zero, 



£^ = 0. 

3.4 Result and a conjecture 

Up to third order in the mass difference we have obtained the ground state energy 



(59) 



E — E —-_ 



h 2 TT 2 



\QmL 2 I mi m 2 



1 + 2— = —B 



L 



h 2 n 2 h 2 n 2 a h 2 n 2 



4mL 2 2mL 2 L 16mL 2 \mi 



m 2 



— -— 1 + 2^ 



(60) 
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to first order in a/L. 

With unequal masses the total mass mi 
connection (l4"5l) between the masses, we have 



m 2 of the molecule is no longer 2m. Using the 



1 



mi + m 2 



1 

2m 



1 



mi + m 2 



1 

mi 
1 



1 

m 2 



777-17772 



4(mi + m 2 ) 



mi + m 2 

Our result (|60|) takes then the simple form 



1 

2m 



1 - 



rn 

77) l 



1 

mi 

777 

m 2 



1 

m 2 



T 



i 

mi 



1 

m 2 



2^2 



-B 



2(mi + m 2 )L 2 



2(mi+m 2 )(L-a) 5 



(61) 
(62) 

(63) 



to first order in a/L. Hence the dead-layer width is a/2, the same value as for equal masses. 

Thus perturbation in the mass difference gives no contribution to the dead-layer width in first, 
second and third order. Although contributions in higher-order perturbation theory cannot be 
excluded a priori, we conjecture that in one dimension the dead layer is determined solely by the 
reduced mass, independent of the mass ratio. 

Intuitively one may perhaps understand a null result in one dimension as follows. For unequal 
masses one particle will have a shorter average distance to the center of mass than the other particle. 
In one dimension the two distances are given by 



(\x 1 -X\) = (\x\) 



m 2 



mi + m 2 



and (\x 2 -X\) = (\x\) 



mi 



mi + m 2 



(64) 



If one believe that the excluded zone is due to the short and the long distance with equal probability, 
the effective width of the dead layer should be determined by the average of the two distances in 
(|64[) . which is independent of the mass ratio. This heuristic argument does not apply in higher 
dimensions. 



4 Concluding remarks 

We have obtained an exact solution for a system of two interacting particles, each of mass m, confined 
to a one-dimensional volume of size L between hard walls. A dead layer is shown to exist in the 
limit L — > oo, and the dead-layer width is precisely equal to the size (|ir 2 — xi\) of the unconfincd 
two-particle system. For finite L the dead-layer width is smaller than this value. We have also 
determined the deformation of the system close to a wall. 

When the two particles have unequal masses, we have determined the dead-layer width by per- 
turbing in the mass difference. Up to third order in the mass difference the dead-layer width is 
unchanged, and we conjecture that in one dimension, and for a given reduced mass, the dead-layer 
width is independent of the mass values. 

Appendix 

The second-order perturbation contribution (|55|) to the energy is 

^ (2) = f f-~-) 2 G, (65) 
8m \777i m 2 / 
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with 



G=(0 



d d 2 

x 







(66) 



dx dX 2 

Now we evaluate this contribution, using the ground state (1261) . 

(X, x\0) = ip(X, x) = A[cosh(u|x| - u ) cos{2vX) - cosh(2wA) cos(w|x| - v)\. (67) 

Since both the operator in (|66p and the wave function (|67[) have even parity in x and in X, we 
may restrict ourselves to positive X and x, and multiply the result by 4: 



G = 4 



dx 



(l-x)/2 



(9 



(9 2 



dA^A,*) — x-^ ij,(X,x) 



A further simplification is that we need the result only to zeroth and first order in the small 
parameter a/L. We use that the parameter u is very large, u ~ L/a to dominating order, and that 
v has the expansion (|36[) . 



1 



a 
I 



a 

L2 



(9 



V' 



The dominating behavior of the normalization constant ([29)) is as follows 



N 2 = \ 



&uv 



L 2 (sinh(2u) + 2u)(sin(2u) + 2v) 
In the wave function (|67[) we use that 



cosh(wa; — u) = cosh(u) [cosh(ua;) — tanh(u) sinh(ux)] ~ cosh(u)[cosh(ua;) — sinh(wx)] ~ |e u e 



(69) 



(70) 



(71) 



since tanh(u) = 1 plus exponentially small terms. The term (|71[) contains the factor e u and it will 
dominate the second term in ip, ea. (|67p . Although all integrations may be performed exactly, it 
suffices for our purpose to approximate the wave function (j67j) as follows 



tP = ±Ne uy cos(2vX), 



(72) 



where y = 1 — x. Thus 



G 



N 2 v 2 / dx e 2u(1 - x) (1 - ux) 



(l-x)/2 



cos 2 (2i;A) dX 



= -N 2 v 2 



? 2uy (l-u + uy) 



y t sinewy) 
4 " 



8c 



dy 



-N 2 v 2 e 2u 



1 

16u 



u 3 sin(2i;) + 3uv 2 sin(2t>) — 2v 3 cos(2w) 



32v(u 2 



2„,2„2u 



-A^t/e 



1 



sin(2u) 



32uv 



using that u = L/a plus exponentially small terms. By means of the expansion (|69[) for v and the 
expression ([70| for the normalization constant A, we obtain to the order of interest 

(73) 







)--£('- 


h2 z) 







By insertion into (|65j) we obtain 
This is eq.(|57P in the main text. 



2^2 



16toL 2 



m 
m i 



m 
m 2 



(74) 
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